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ABSTRACT 

We discuss the probability distributions of Fourier modes of cosmological den- 
sity fields using the central limit theorem as it applies to weighted integrals of 
random fields. It is shown that if the cosmological principle holds in a certain 
sense, i.e., the universe approaches homogeneity and isotropy sufficiently rapidly 
on very large scales, the one-point distribution of each Fourier mode of the density 
field is Gaussian whether or not the density field itself is Gaussian. Therefore, 
one-point distributions of the power spectrum obtained from observational data 
or from simulations are not a good test of whether the density field is Gaussian. 

PACS number(s): 98.80_k, 98.80.BP, 98.60.Eg 
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I. INTRODUCTION 



In the standard cold dark matter (CDM) model, primordial fluctuations are 
generated during inflationary stage, and the simple linear perturbation theory pre- 
dicts Gaussian fluctuations [1]. On the other hand, some alternative inflationary 
theories [2] [3] and topological defect mechanisms [4] [5] involve nonlinear pro- 
cesses which generate non-Gaussian perturbations. Thus it is interesting to test 
the consistency of a theory with observations in this regard. Besides analyzing the 
spatial distribution of galaxies [6] [7] and the pattern of microwave background 
radiation anisotropies [8] , some research groups considered the one-point distribu- 
tion of the power spectrum of the density field as a test of the Gaussianity [9] [10] 
[11] 12] [13]. It is found that for the QDOT IRAS redshift survey, the distribution 
is almost exactly an exponential distribution which is consistent with the Gaus- 
sian distribution for both the real and imaginary parts of a Fourier mode [13]. 
Suginohara and Suto [14] did simulations up to very nonlinear stages, analyzed 
the one-point distribution for both the power spectrum and the phase of a Fourier 
mode, and also got results consistent with Gaussian distributions. On the other 
hand, we know that in highly nonlinear stages the one-point density distribution 
in real space is very non-Gaussian. This raises the question, namely, what can 
one really conclude from the statistics of the one-point distribution of the power 
spectrum? 

The Fourier transform of the cosmological density field is basically a weighted 
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sum of the densities over the space. If the universe consists of a large number 
of dense clumps (e.g., galaxies, clusters), and those clumps are independent of 
each other, then by virtue of the central limit theorem, the real and imaginary 
components of each individual Fourier mode are Gaussian distributed although 
the density distribution itself could be highly non-Gaussian. Kaiser and Peacock 
[10] generated some test catalogs, and analyzed the distributions of Fourier modes. 
They found that even if only a few independent clumps existed in a test catalog, 
the one-point probabilities of individual Fourier modes would be well described 
by Gaussian distributions. 

On the other hand, if there are some correlations on scales larger than the 
typical size of a clump (as most of the cosmological models predict), what is the 
one-point probability of an individual Fourier mode? In this paper, it is shown 
that if the random density field satisfies certain mixing conditions, the central 
limit theorem guarantees that the one-point distribution of an individual Fourier 
component is Gaussian whether or not the field itself is Gaussian. 

The paper is organized as following: In section II, we introduce some concepts 
and discuss an application of the central limit theorem to weighted integrals of a 
random field. In section III, we consider the Fourier transform of a random field 
and see how it satisfies the conditions in section II. Section IV contains further 
discussion and conclusions. The appendix provides an approximate estimate of 
the mixing rate, which is introduced and used in the section II, in terms of the 



two-point joint probability. 
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II. THE CENTRAL LIMIT THEOREM 



In this section, we discuss a central limit theorem for weighted integrals of 
random scalar fields in terms of constraints on the field moments, correlations 
and the weighting functions. The theorem and its context is clearly discussed by 
Ivanov and Leonenko [16]*. 

Let us first consider a random variable 5. Without losing generality, we assume 
the average < 5 >= and 5 G [ft, /„], where / ; G -R 1 and /„ G R 1 . Let p(S) be 
the probability density function of 5 : and a set A s be {5 : f t < a l < 5 < b l < /„}. 
Then the probability that 5 taking values on the set A s is 

P(A 5 )= [ bl p(5)d5 , (2.1) 

and the normalization condition is 

[ U p(5)d5 = l . (2.2) 
Jfi 

Now consider a random scalar field 5(x) satisfying < 5(x) >= 0, x G -R 3 . We 
assume that the probability density function and the value range of 8{x) are the 
same for all x G -R 3 . Let A x and A 2 be two regions in R 3 . The distance between 
A x and A 2 is defined as 

r = min|x — y\, x G A 1 , y G A 2 . 

In order to apply the central limit theorem, we need to introduce the concept 

of the mixing rate [17]. Let us denote positions on A l5 A 2 by {x„i = 1,2, ...}, 
* See the reference 16, Chapter I. 
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{yjij = 1, 2, ...}, respectively. Let A be a set generated by the random variables 
{5{xi), i = 1, 2, ...} and 

a* < S(xi) < b u a, > f u b, < f u , i = 1, 2, ... , (2.3) 
and let B be generated by the random variables {S(yj),j = 1, 2, ...} with the form 

Cj<S%)< ej , Cj >f h e,<f u , j = 1,2,... . (2.4) 
Then the probability P(A), P{B) and P(AB) are 

P(A) = {II, T ^(f,)} pWi), 8{x 2 ), ...} , (2.5) 

J CLi 

P(B) = {II, P <W(j&)} p[5(yi), •••] , (2-6) 

and 

p(ab) = {n iJi=li2 ,... r^(^) r^)}^^,^),...,^!),^),...] . 

(2.7) 

respectively, where p[5(xi), 8(x 2 ), ...] is the joint probability density function of 
{5(xj),z = 1,2,...}, p[S(yi),S(y 2 ), ■■■] is the joint probability density function of 
{8{y ),j = 1,2,...}, and p[8(x i), 8(x 2 ), 8(y 2 ), ...] is the joint probability 
density function of {5(5;,), £(%•), i = 1,2, = 1,2,...}. The dependence of 
the random variables defined on A x and A 2 is characterized by the Rosenblatt 
dependence rate [17] which is defined as 

a{A u A 2 ) = max \P(AB) - P(A)P(B) \ , (2.8) 
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where the maximum is taken over all possible values of {a i7 6, G [/;, /„], c,, e,- G 
[/«) /«]) * = 1) 2, j = 1, 2, ...}. The mixing rate a(r ) is defined as 

a(r Q ) = maxa(A 1; A 2 ) , (2.9) 

where the maximum is taken over all pairs of A 1 and A 2 lying at a distance at 
least r . It is clear that in general the mixing rate a(r ) is not only related to the 
two-point but also to all higher order correlations. 

With these definitions, we now discuss the central limit theorem, and introduce 
a column vector s v (9) = [sy\9), Sy\0)]\ where t stands for the transpose, with 
components defined by 

s$0)= [ g { v\x,6)8(x)dx , (2.10) 
Jv 

where g v (x,9) = [gy\x, 9), ...,gy\x, v)] 1 is a column vector of functions, and 9 
G 0. In the situation we are going to discuss in the next section, 6 = R 3 . 

We assume that the following conditions are satisfied: 

(i) . As V — > oo, there exists a function cr 2 (9), taking on values in a set of 
positive definite r x r matrices such that 

lim [< s v {9)s v {9) > -a 2 (9)] = (2.11) 

— * 

uniformly in 9 G 0. 

(ii) . For any i G {1, ...,r}, the function gy\x, 9) satisfies the condition: 

For V — > oo, there exists a constant > (possibly depending on i) such 
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that 

max \g^(x,6)\<k^/W . (2.12) 
eee.xev 

(iii) . There exists a constant 7 > such that 

< \8{x)\ 2+1 >< kx < 00, xeR 3 . (2.13) 

(iv) . The mixing rate a(r ) of the random field 8(x) satisfies 

«(r ) < /e 2 r - 3 " £ , e 7 > 6 , (2.14) 

where k 2 is a constant. 

— * 

Then as V — > 00, the random vector s y (6>) has, uniformly in 0, an asymp- 
totically multidimensional normal distribution N r (0, c 2 (9)) with zero mean vector 
and covariance matrix a 2 (9)*. 

In the next section, we will discuss the Fourier transform of a cosmological 
density field using this central limit theorem. 



* See the reference 16, p41, Theorem 1.7.5 
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III. FOURIER TRANSFORM OF A RANDOM FIELD 

The Fourier transform of a random field 8(x) in a 3-dimensional space is 
defined as 

8(k) — lim — = [ expOik ■ x)8(x)d 3 x . (3-1) 
v^co yjy Jv 

Here we assume 8{x) is a homogeneous, isotropic and real random field, which 
satisfies < 8{x) >= 0. 

The Fourier transform 8{k) can be written as a sum of the real part A(k) and 
the imaginary part iB(k), where A(k) and B(k) are defined as 



and 



A(k) = lim Ayik) , (3.2) 



B(k) = lim B v (k) , (3.3) 



respectively, where 



A v = —= [ cos(k ■ x)8(x) 
W Jv 



and 



By = —== I sin(/c ■ x)8(x) 
\/V Jv 



— * — * — * 

We consider the vector s v (k) = [A v (k), B v (k)]. In terms of the notations of 
section II, we have r = 2 and 



g ( v 1 \x,k) = ^=cos(k-x) , (3.4) 
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9 ( v (x, k) = -j= sin(&; • x) . (3.5) 

We discuss the conditions of the central limit theorem in the section II one by 
one. 

(i). To evaluate the covariance matrix a 2 (k), three elements should be con- 
sidered: 

alAk) = lim < A v (k)A v (k) > , 
aUk) = lim < B v (k)B v (k) > , 

and 

al 2 (k) = hm < A v (k)B v (k) > . 

First let us consider o\ x . 
a ii(k) = lim / / d 3 x 1 d 3 x 2 — 7= cosCk-xA — ]= cos(k-x 2 ) < 5(x 1 )5(x 2 ) > ■ (3.6) 

V-KX> Jy J V y/y yfy 

Because 5(x) is homogeneous, isotropic and < >= 0, < S(x 2 ) >= 0, the 

correlation < 5(xi)5(x 2 ) >= £(|ici — x 2 \) only depends on \x 1 — x 2 \. One can 
readily show that 



a li(k) = Jlim f d 3 x cos(k ■ x)£(\x\) . 

2 V^oo Jy 



We know that the power spectrum averaged over the ensemble of realizations 
is related to the correlation function £(x) by 



< P(k) >= lim / d 3 xexp(ik ■ x)^(x) 

V^co Jy 
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Since S(x) is homogeneous, isotropic and real, it follows that 
< p(k) >=< P(-k) >=< P(k) > and f (£) = £(\x |). Thus 

lim / d 3 xcos(k ■ x)£(\x\) =< P(k) > , 



and 



lim / d 3 x sin(/c • x)£(|a;|) = 



Then 



^(fe) = i<P(fe)> • (3.7) 
Similarly one can show that 

a 2 22 (k) = l<P(k)> , (3.8) 

and 

a 2 12 (k)=0 . (3.9) 

For any reasonable cosmo logical density field, the average power spectrum 

< P(k) > exists. Therefore the condition for the existence of the covariance 
matrix cr 2 (k) is satisfied. 

(ii). From (3.4), (3.5), | cos(/c • x)\ < 1, and | sin(/c • x)\ < 1, we can always 
choose k 3 = 1 so that condition (ii) in section II is satisfied. 

Condition (iii) of section II constrains the statistical properties of the field 
5(x). The existence of moments higher than the second is not a problem for most 
models of cosmological density field. 
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Condition (iv) of section II concerns the behavior of the mixing rate a(r ). In 
general this mixing rate is related to all order of correlations and thus is difficult to 
be calculated precisely. Under certain approximations, however, we can estimate 
the behavior of a(r Q ). In the Appendix, we make such an approximation: Between 
two well separated regions A x and A 2 , only two-point joint probabilities are con- 
sidered, while within each region full correlations are taken into account. We 
argue in the Appendix that under plausible assumptions, the reduced two-point 
joint probability of Si at x t and Sj at pj can be written as 

oo 

m,6j,\2<-dj\)= E M^S^FTM-yjl), (3.10) 

n=no 

where F i:j (\xi — y*j\) <C 1 when \x { — y 3 \ is very large. A n (8 u Sj) is independent of 
Fij(\xi — and n is the power index such that A n vanishes for n < n Q . Here 
n a > 1. 

Under these approximations 

a(r )~C' max \F?;(\x t - (3.11) 

where C' is a finite constant, and the maximum is taken over all pairs of (xi, yj) 
with \xi — yj\ > r . The condition (iv) of section II requires that 

«(r ) < k 2 r- 3 ~ e , e 7 > 6 . (3.12) 

Then in order the inequality (3.12) to be satisfied, we need 

rnax \F^(\Si - yj\)\ < C 14 r" 3 - , (3.13) 
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where C u is a finite constant. As discussed in the Appendix, for a non-Gaussian 
field which is a local functional of a Gaussian field, F io {\xi — y \) is the two- 
point correlation function of the Gaussian field. The power index n in (3.10) 
depends on the specific form of the non-Gaussian field. For example, for a field 
$ = (p 2 — < (p 2 >, where is a Gaussian field and < (p 2 > is the average of 4> 2 , 
the power index n is 2. If the two-point correlation function of a Gaussian field 
satisfies condition (3.13) (in this case n = 1), all non-Gaussian fields which are 
local functionals of the Gaussian field satisfy the mixing rate condition. 

Reduced higher order joint probabilities are related to higher order reduced 
correlations. Observations indicate that higher order of reduced correlations drop 
faster than the two-point correlation function [15] [18] [19]. Thus it is a reasonable 
approximation for a cosmological density field to ignore reduced higher order joint 
probabilities. 

As a case where the mixing rate condition is violated, we consider a random 
field S(x) which is a convolution of two Gaussian fields (fi 1 and 4> 2 

5(x) = J d'xMx^ix -xj) . (3.14) 

Clearly, for each x G -R 3 , S(x) is related to the fields and 4> 2 in the whole space. 
Thus the probabilities of the field S(x) is always related to the joint probabilities 
of all and all (fi 2 in the whole space. It is not surprising that some complicated 
correlations for the random field 5(x) exist and the mixing rate condition is not 
satisfied. In fact, we know that the Fourier transform of S(x) is the product of 



14 



the Fourier transform of </>i and </> 2 



8(k) = , 



(3.15). 



— * 

so 8(k) is not Gaussian distributed. 



If a field S(x) satisfies the conditions of the central limit theorem, then the the 
joint distribution of the real and imaginary parts of its Fourier transform, namely, 
A(k) and B(k), is simply 



p(A, B)dAdB 



n < P(k) > 



exp 



A 2 + B 2 
< Pjk) > 



dAdB 



(3.16) 



If we change the variables from A(k) and B(k) to P{k) and 6>£, where 9^ is 

— * 

the phase of the k mode, we get 



P 



P(k),e s 



— exp 



P(k) 



< P(k) > 



dPddt 



k -i 



(3.17) 



that is, P(k) is exponentially distributed, and 9$ is uniformly distributed. 
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IV. DISCUSSION 



In this paper we showed that as long as the mixing rate of the density field 
drops fast enough at large distances, the central limit theorem guarantees normal 
distributions for the real and imaginary parts of a Fourier mode of the field. 

Reduced higher order joint probabilities are related to the reduced higher 
order correlations. Observational evidences show that on large scales, these cor- 
relations do drop faster than the two-point correlation function does [15] [18] [19]. 
Therefore for the mixing rate, it is reasonable to only consider the two-point joint 
probabilites between two well separated regions. Under such an approximation, 
we estimate the mixing rate. If the reduced two-point joint probability of S { at x { 
and Sj at y 3 depends on a function F^Qxi — y 3 \) in such a way 

oo 

p(5 t ,5 v \x t - y 3 ]) = A n{&i^j) F ?j{\xi ~ Vj\)y 

n=no 

where \F^(\xi — y 3 \)\ <C 1 when \x t — y 3 \ is very large, and A n {5 u S 3 ) is independent 
of 

~ Vj\)i th e mixing rate has the same behavior as 
^ max 1^7(^-^1)1, 
where C is a finite constant. Thus if 

max 1^(1^- yj\)\ <C 14 r " 3 - 

\xi-Vj\>r 

at large r , where e > 0, the mixing rate condition is satisfied. The condition on 
max|^_^.|> ro \F£°(\xi — y 3 \)\ is equivalent to requiring that 

max p(8i, 8 3 , \xi - y 3 \) < C 15 r 3 ~ £ . 

\Xi-Vj\>r 
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This is a sensible requirement, and is easily satisfied by a cosmo logical density 
field. 

Statistical properties of the large-scale structure in the universe provide im- 
portant clues for understanding the universe. Theoretically, to analyze statistical 
properties of a random field, one must generate many different realizations of the 
field and calculate the statistical quantities over these realizations. Practically, 
since there is only one universe, we have to use the " ergodic" theorem in the 
position space to deal with cosmological density field. Therefore it is important 
to have what is in some sense " a fair sample " . The basic requirement for "a fair 
sample" is that quantities within the sample should represent the quantities of the 
whole universe without being strongly biased by the correlations of the quantities 
within the sample. If the universe is statistically homogeneous and isotropic at 
very large scales, all correlations within "a fair sample" should have no signifi- 
cant effect on the statistical quantities and can be ignored. For such samples, the 
Fourier modes of the density field are Gaussian distributed whether or not the 
density field itself is Gaussian. Since in reality one can never be certain whether 
or not we really have a fair sample, we must be cautious to interpret the statistical 
results on the distributions of the Fourier modes. If non-Gaussian distributions of 
individual Fourier modes are found for a set of samples and this non-Gaussianity 
disappears when larger samples are taken, it is most likely that the original sample 
is not a fair sample, rather than that there is a genuine non-Gaussian behavior. 
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If the non-Gaussianity persists when we go to larger and larger samples, then one 
might begin to believe that there is some physics in it. 

One class of non-Gaussian models are topological defect seeded models, e.g., 
cosmic strings [5], or textures [4]. In those models, the density perturbations are 
due to randomly distributed seeds, and each seed has a density profile. If the 
density profiles of all seeds are the same, and are denoted by f(x — where x { 
is the position of a seed, the density field can be written as 

N 

S(x) = J2f(x-x i ), (4.1) 

i=l 

where the sum is over all the seeds. We can write (4.1) as 

S(x) = J dx'g{x')f{x-x'), (4.2) 

where 

N 

g(x')=J25^(x'-x t ) 7 (4.3) 

»=i 

and 5 (3 \x' — Xj) is the three-dimensional Dirac 5-function. Then the Fourier 
transform of S(x) is 

F 5 (k) = F g (k)F f (k), (4.4) 

where F g (k) is the Fourier transform of the function g, and F f (k) is the Fourier 
transform /. Since / is not a random function, F f {k) is not random. The distri- 
bution of F s (k) depends on the statistics of F g (k). We have 
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If the seeds are uncorrelated and N — > oo, F g (k) is Gaussian distributed and so 

— * 

is F s (k) by the central limit theorem for the discrete summation[20]. Certainly, 
if within an observational region only one or two seeds exist, one may get some 
non-Gaussian signal in the distributions of the Fourier transform [20]. But in 
that case, one needs to go to larger samples in order to get statistically solid 
conclusions. 

In principle, the appropriate way to analyze the distributions of the Fourier 
modes (power spectrum) is to generate a large number of realizations of the cos- 
mological density field. For each realization, the Fourier transform is made, and 
the average power spectrum < P(k) > for each k is obtained by averaging over 
all realizations. Then the distribution of P(k)/ < P(k) > for each k is analyzed. 
Practically, since we only have one universe, and can do limited number of sim- 
ulations, some kind of "ergodic" assumptions are always used when analyzing 
observational or simulation data. For example, in [12] and [13] the distributions 
of the power spectrum in /c-space are analyzed, i.e., to cumulate the numbers of k 
which have the same value of Pj < P >, where < P > is the average power based 
on some physical considerations. Then the distributions of those numbers against 
P/ < P > are analyzed. In these analyses, we should be very careful to use an 
appropriate < P >. For example, as pointed out in [13], because of the small 
scale clustering, the average < P(k) > is larger than the purely random walk esti- 
mation for large k. If one uses < P > estimated by ignoring the clustering effect, 
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one may get a power spectrum distribution which is different from the exponential 
distribution expected from Gaussian distributions of Fourier modes. This differ- 
ence is caused by using an inappropriate average power spectrum, rather than the 
existence of some genuine non-Gaussianity. 

By the above discussion, we see that the one-point distributions of Fourier 

modes are not good statistical tests of the Gaussianity of the density field. One 

possible way to distinguish Gaussian and non-Gaussian statistics is to go to high 

order correlations. For example, we can consider the two-point correlation of 

the power spectra F(k,k') =< P(k)P(k') >, which is related to the four-point 

correlation of Fourier components [13]. For a Gaussian random field (j)(x), 
F(k,k') = < P(k) >< P(k') > 

+ 5 {3) (k + k') <P(k)><P(k)> (4.1) 

+ 5 {3) (k - %') < P(k) >< P(k) > 
If k + k' ^ and k - k' ^ 0, then 

F(k, k') =< P{k) >< P(jfe') > (4.2) 

i.e., for a Gaussian field, the correlation of the power spectrum for different k- 
mode is simply the product of the power spectrum of each mode. 

For a non-Gaussian field, it is expected that there are some additional terms 
in the correlation of the power spectrum. To illustrate this point, as an example, 
we consider a field ((x) which is the square of a Gaussian random field 4>{x). 

m = <m (4.3) 
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Detailed calculations show that if k ^ —k' and k ^ k', 
F{k, k') = < P{k) >< P(jfe') > 

+ s/ 1^ < \tfymk - feon^i + k')\ 2 \<t>(k- k' - h)\ 2 > 

+ 8 j^y 3 < \<P(h)\ 2 \<P(k - fcon^fex - £')I 2 I<K£+ # - > 

(4.4) 

It can be seen that there are some extra terms added up to < P(k) >< P(k') >. 
Generally these extra terms depend on the specific form of the non-Gaussian field. 
Therefore at this level, it is possible to distinguish Gaussian and non-Gaussian 
fields. 
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Appendix 

Estimation of the Mixing Rate 

In this appendix, we discuss the approximate behavior of the mixing rate. 

Consider \P(AB) — P(A)P(B)\ as discussed in section II. If the random vari- 
ables in Ai are completely independent of the random variables in A 2 , then 
P{AB) = P(A)P(B), and \P{AB) - P(A)P(B)\ = as expected. If the two 
sets of random variables are weakly correlated, we could only consider the first 
order corrections to P{AB) in addition to P(A)P(B). We consider a statistically 
homogeneous and isotropic random field. For two random variables Si at x t and 
Sj at yj, their joint probability is p(Si, Sj, \x t — yA). Then 

pi&uSj^Xi-y^) = p(S h S jl \x i -y j \) -p{5 i )p{8 j ) (AI) 

is the reduced two-point joint probability which represents the contribution to 
p(Si,Sj, \xi — yA) due to the correlation between S l and Sj. The first order ap- 
proximation of P(AB) — P(A)P(B) means that between two well separated re- 
gions, only reduced two-point joint probabilities are considered. Let us consider 
the contribution from the correlation between Si G A : and S 3 E A 2 . After 
5, 6 Aj having been picked out, since Si is correlated with the rest of the vari- 
ables within Ai, the joint probability for the rest of the random variables within 
A 1 should be the conditional joint probability with Si fixed, which is denoted by 
p(S(xi), S(x i+ i), ...\S(xi)). Similarly, the joint probability for the rest of 
the variables within A 2 should be p(S(yi), 5(^ 7 _ 1 ), S(y j+1 ), ...\S(yj)), the condi- 



(A2) 
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tional probability with 5j fixed. Therefore the contribution to P(AB) — P(A)P(B) 
from the correlation between 6^ and Sj has the form 

p ij={ I dS(xi) [ dStij^pfaSj^Xi- yj\) 

f hpi -> [ e "i -> 

^■p i =i,...,i-i,i+i,...,g j =i,...j-i,j+i,... / dS(x p .) / dS(y q .) 

Ja p . Jc q . 

p(<y(fi), ...|<y(fi)) 

The first order correction to P{AB) — P(A)P(B) should include contributions 
from all pairs of random variables between A : and A 2 . When r is sufficiently 
large, approximately we have 

P(AB) - P(A)P(B) ~ £ P ^ ( A3 ) 

i,j=l,2,... 

where P y has the form (A. 2) 

Since we are interested in the behavior of the mixing rate with large separation 
between the two regions, we would like to separate the correlations between Si and 
Sj, which depends on y*j\ for a statistically homogeneous and isotropic random 
field, from the statistics for each of the random variables. First, let us consider a 
non-Gaussian field S(x) which is a local functional of a Gaussian field <j>(x), i.e., 
5(x) = 5[4>(xj\. The joint probability between 8(xi) and 8(yj) is related to the 
joint probability of 0(a?j) and 4>{Vj)- We assume that for a fixed value of S(x), 
there are m values of (ft k \x), (k = 1, m) corresponding to it, it then follows 

m 

p[S(x i ),S(y j )} = E P[^M^\m 
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where the factor {l/|£/5[0(x i )]/d[0(x i )]|fe i }{l/|d5[0(y :7 -)]/£i[0(y J -)]|fe i } comes from 
the Jacobi transformation taking value at (fei) (^) and ft ki \yj). It is well known 
that the two-point joint probability p[<j>{xi),<j>(y*j)] of a Gaussian field 4> with 
< </> 2 >= 1 can be expanded as 

P [mum] = E ^v^^ ( " +i) [0(^)]^ (rt+1) [0(%-)], (^-5) 

where — y,-|) is the two-point correlation function of the Gaussian field 4>(x), 

and 

$ ( " +1) [<^)] = ^exp[-^(fO/2]*W(2,)], (^-6) 
where H n is the Chebyshev-Hermite polynomials *. Then (A. 4) can be written as 

oo 

p[S(x i ),S(y j )] = J2C^i-yMn[S(x i ),S(y j )}, (A7) 

n=0 

where 

1 m r 1 

+ ( - )& n+1) \6 {k > ) (v)}\ 

+ { \dsm)]/dm)\w [ [Vj)l ) 

The reduced two-point probability is 

oo 

m*i),S%)] = J2Q(\^-yj\) A n[S(x i ),S(y j )}. (A9) 

71=1 

For a specific non-Gaussian field, the index n of the first nonvanishing A n [S(xi), 
can be greater than 1. We can write (A. 9) as 

oo 

p[6(2MPj)] = E Qdx.-yjDAM^^m, (A.10) 



{A.i 



See the reference 16, Page 55 
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where n is the index such that all A n (n < n ) vanish. The dependence of 
p[S(xi), S(yj)] on \x t — yj\ is thus shown explicitly in (A. 10). 

Based on the above considerations, we make several assumptions: 

(i) . p(Si, Sj, \xi — yj\) depends on a function F i3 (\xi — yj\) and can be written as 

oo 

m^ r \^-y,\)= e M^^FrM-yji), (ah) 

n=no 

where A n (S l ^ Sj) is independent of F i3 (\xi — y 3 |). When F y (| Xi yj 

I) = 0, (A.ll) 

implies that 5^ and Sj are independent of each other and p(d i ,d J ,\x l — %|) = 
0. For a non-Gaussian field which is a local functional of a Gaussian field 0, 
Fyda^ — yj|) = — %|) as shown above. The random variables can have 

value ranges extending to infinity, so it is necessary to keep all the terms in the 
expansion (A.ll). 

(ii) . When \x t - y 3 \ is large, F i3 (\x, - y \) < 1. 
By (A.ll), Pij in (A. 2) can be written as 

OO 

Pij= E W-ftK (A.12) 

n— no 

where 

C ntj ={ / * i) 1 ' d5(y 3 )A n (5 t , 5,) 

I J (Xi J Cj 

f b pi f e 1j 

n Pi=1 ,... ;i _ M+1 ,.. / d5(x Pi ) / cM(y g .) 

(A13) 

p(<y(fi), <y(£i_i), <y(x i+ i), ...|<y(fi)) 
x p[*(i7i), 5(%-0, <y(y j+ i), ...|<y(?7i)] k 
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is a function of (x u ...,a u b u ...,a t ,b t , ...) and (y u ...,c 1 ,e 1 , ...,c 4 , e t , ...). 

Since C„ i3 is finite, by the assumption (ii), we only consider — y 3 \) and 

ignore all the other higher order of F ij (\x l — y*j\). Then 

P(AB) - P(A)P(B) ~ (I** - &|)C n0y (A14) 

For a specific random field, F tj may have a certain value range. It is clear that 
Fij = must be an allowable value. 



We also have the following relation: 

\P{AB) - P(A)P(B)\ ~ | 52F?{\2 t - y 3 \)C n ^ \ < £ \F^(\x t - yj\)C no J, 



El^ (l^-%l)^oJ < .max -%)HEl^oJ, (A15) 

where the maximum is taken over all possible Xi G A : and G A 2 , and 
max |P(AB) - P(A)P(B)\ ~ max | £i^(|f< - fcDO^ I 

(A16) 

< .max |^7(|^ -y,|)|max(Eia J), 

1,3 

where the maximum on the left hand is taken over all possible integration ranges. 
On the right hand of (A. 16), the first maximum is taken over all possible x l G A x 
and y*j G A 2 , and the second maximum is in the same sense as the maximum on 
the left hand. 

We know that \P(AB) - P(A)P(B) \ < 2 regardless of both the specific values 
of (within the value range) and the specific integration ranges for the random 
variables. Then the sum J^i.j \C n0ij I should converge. Concerning the approximate 
behavior when the separation between Ai and A 2 is sufficiently large, we have 

max \P(AB) - P(A)P(B)\ ~ C _ max \F^{% - (A 17) 
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where C is a finite constant. By the definition of the mixing rate a(r ) in (2.9), 
we have 

a(r )~C max 1^(1*, - &|)|, (A.18) 

\Xi-yj\>ro 

where C' is a finite constant and the maximum is taken over all pairs of y } ) 
with \xi — Vj\ > r . 



